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Abstract 

We compute CP asymmetries of one-particle inclusive B® — ► D S X decay 
rates. We find that the weak angle 7 could in principle be extracted from 
these asymmetries although they are very small. Taking the decay width 
difference into account, we find that the time integrated CP asymmetries 
are of the order A = 1.4 10 -4 sin(7). Some large uncertainties remain due 
to the lack of experimental data. 
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1 Introduction 



Some time ago it was proposed by Beneke, Buchalla and Dunietz |J to study 
mixing induced CP violation in B mesons by using asymmetries of inclusive de- 
cay rates. This method although theoretically clean suffers from an experimental 
difficulty since it supposes a completely inclusive measurement. In order to solve 
this problem, we propose to study CP asymmetries of one-particle inclusive de- 
cays. In this paper, we shall consider CP asymmetries in the B® system. This 
can also be done in the Bd system 0. We shall study time integrated and time 
dependent CP asymmetries of one-particle inclusive B° s — > D S X decays. The 
technique involved is nevertheless theoretically not as clean as the one used in 
[0 for inclusive decays. Nevertheless, the measurement of CP asymmetries of 
one-particle inclusive decay widths would be clean. 

We shall not neglect the decay width difference in the B° s system, since it 
is sizable although recent calculations || show that the decay width difference 
in that system could be smaller than previously expected. A method |4]] was re- 
cently proposed to calculate the decay rates appearing in these asymmetries. The 
predicted decay rates are compatible with current experimental knowledge and 
we can therefore be confident that the results we obtain for the CP asymmetries 
will be meaningful. 

We shall concentrate on asymmetries involving one-particle inclusive decays of 
a B° s meson. We start by introducing our notations, in section |^ we shall develop 
the formalism needed for time dependent CP asymmetries in one-particle CP 
asymmetries. We then present and discuss our results and conclude. 



2 Definitions 



In this section we shall introduce our notations. We basically use notations similar 
to those introduced in . The proper time evolution of an initial pure B® or B® 
reads 



B° sphys {t)) = g+(t) | B° s ) - lg_(t) | B° s ) 

(t)) = -^g-(t)\B s )+g4t)\B° s }, 



(1) 



dO 

s phys 



where the time dependent functions 

-iMt-lrt 



g+(t) 



, Art AMI , Art AMI 

cosh cos h i smh sin 



(2) 



and 



9-® 



iMt~\Yt 



, Art AMt , Art AMt] 
smh cos h i cosh sin (3) 
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describe the particle anti-particle mixing. The mass difference AM and the width 
difference Ar between the neutral B° s mesons are given by 

AM = M H -M L (4) 

Ar = r H -r L . 

The off diagonal term of the mass- mixing matrix is given by M\i + iV\2- We 
define the quantity 

p 2{M 12 -i/2T 12 ) |M 12 | V 2 J' M 12 v; 

The second expression for q/p in equation (|5|) is valid to first order in the small 
quantity T 12 /M 12 = 0{ml/m 2 t ). 



3 CP asymmetries in one-particle inclusive de- 
cays 

In this section, we shall consider CP asymmetries of rates of one-particle inclusive 
— > D S X decays. The time dependent CP asymmetry is defined by 

A m r(B°(t) -> D S X) - T(B° a (t) -> D S X) 

•* CP{ ) r(B°(t) -> zvO + r(£°(t) -> d s x) " 1 J 

In the following, we shall neglect all effects due to direct CP violation. The time 
dependent decay widths therefore read 

T(B° s (t) - D S X) = \g + (t)\ 2 r(B s -> D S X) + | V(t)| 2 r(£ s ° -> D S X) (7) 

-2Re (^(t) 

and 

T(B° s (t) - D S X) = \g + (t)\ 2 T(B° a - D S X) + | V(t)| 2 r(5 s ° - (8) 

-2Re (gl{t)l 9 -(t) T V «) . 

The decay widths appearing in these formulas were calculated in It remains 
to compute the transition matrix elements of AB® = 2 given by 

Tg ^ = 1^1 J d4>o s E(27T) 4 5 4 (P B o - P Ds - P X ) (9) 



2m B o * ., x 
(B° s \H eff (x)\D s X)(D s X\Hl ff (0)\B° s ) 



and 

T*: B ° S = ^J d'x J d<j> Ds £(2vr) V(P B o - P Da - P X ) (10) 

s X 

(B° s | H eff (x) | D S X) (D S X | R\ n (0) | B° s ) , 

where dcf)£, s is the phase space of the D s meson. The part of the effective Hamil- 
tonian which is relevant for this work is given by 

H eff = ^V; b V cs (C x (bc) v _ A {cs) v _ A + C 2 (bT a c) v „ A (cT a s) v _ A ) (11) 

+^ V aV us (C x (bc) v „ A {us) V - A + C 2 (bT a c) V - A (uT a s) V - A ) 

+^ v : b V cs {Ci{U) v ^ A {cs) v ^ A + C 2 (bT a u) v „ A (cT a s) v _ A ) 

where (qiq 2 )v-A stands for (gi7 M (l — 75)92), the T a matrices are the SU(3)c Gell- 
Mann matrices and C\ and C 2 are the Wilson coefficients. The matrix elements 
T D 3 s " and T^ s 3 can be parametrized in the same way as the wrong charm 
decay in ||. Applying the Fierz transformation we can rewrite the operators into 
the following form (bs)v - A {cc)y _ A , {bs)v- A {uc)y_ A and {bs)y_ A {cu)v- A - Using 

B° B° 

the large Nq limit, we can f actor ize T h s s 

Tg m (M 2 ) = £ ^-^C^VSY^V^ (12) 
J d A QK^( PBh M,Q) J d4>D,P${Q) 

with 

K^(p B o, M, Q) = E(2vr) 4 5 4 (p B o -p x -Q) (13) 

x 

(B° s (p B o)\(h,(l - l5 )s)\X)(X\(b lu (l - l5 )s)\B° s (p B o)) 

and 

P$(Q) = E(2^ 4 (Q - Pd s ~ Px>) (14) 

X' 

(01(^(1 - l5 )q)\D s (p Ds )X')(D s (p Ds )X'\(q lu (l - l5 )c)\0). 

We get a similar expression for ^ a (M 2 ). The nature of the quark q depends 
on the operator under consideration. The tensor K^^pso, M, Q) is fully inclusive 
and can be parametrized using the decay constant f B o of the B° s meson and 
Pjw (Q) * s ver Y similar to the expression we had encountered in the calculation of 
the wrong charm decay width [TJJ. It involves a projection on a state containing a 



D s meson. A priori we do not know how to calculate that kind of matrix element, 
a solution is to contract the spinor indices as in the parton calculation and to 
multiply this tensor by a channel dependent form factor f q which has to be fitted. 
We obtain 

P^(PD,Q) = 2n5((Q-p D ) 2 -m 2 q )Ti{^ Dlll (Q-^ D ) lu }f q . (15) 
We can reproduce the inclusive case by setting f q — 1. 



3.1 CP asymmetry T(B° s (t) D+X) vs. T{B° s {t) -> DjX) 
The CP asymmetry reads 



A(t) 



r(gg(g -> DfX) - T(B° a (t) -» 
r(fl°(f) - D+X) + T(B°(t) - D-X) • 



(16) 



We have two quark transitions contributing to T + s , namely b — » ccs interfering 

U s 

with itself and b — > ucs interfering with b — > cus. These matrix elements can be 
modeled in the same way as it was proposed in |4j] . It corresponds to a rescaling 
of the parton calculation. The complication due to isospin symmetry which led 
to the introduction of strong phases in || is not present in the system. We 
can use the same parameterization for the AB® = 2 as the one we have used for 
the wrong charm decay widths [Q. The parton calculation was performed in M. 
We obtain 



T 



D+ 



Gl '" 1 f B oM*v cb v: s v ub v: 



24nM B o 

Cftl-z)B + C'i(l -f 2:) 

G%m 2 f B o 



24tcM b o 



-M 2 B0 {v c * s v cb ) 2 VT^ 



(1 - 4z)C(B + (1 + 1z)C{-^-2B l 



mi 



J 7 , 



in the leading order of the short distance expansion and where T and Q are decay 
channel dependent non-perturbative form factors. We shall define JF and Q later. 
The parameters B and B$ are the bag factors and z is defined by 



m. 



m 



2 ' 



We have neglected the penguin operators. We also have two operators contribut- 
ing to T ' 3 , as in the previous case b — > ccs interfering with itself and b — > cus 
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interfering with b — > wcs. In the leading order of the short distance expansion we 
then have 



T 



D7 



24nM B o 

2 



/ioM^K5K s K*6K s (i-^) 



C?(l - z)B + C?(l + 2z)^-2B s \Q 

2 ™2 /2 



Gpfn b f B o 
24nM B o 



Mlo(v; b v cs ) 2 VT^ 



1 - Az)CfB + (1 + 2z)C?— 



m 



(19) 



where JF and Q are decay channel dependent non-perturbative form factors. The 
function T was given in M. We had 



b°d: 



4/, 



(20) 



where / = 0.121 was fitted in 0. The value of this parameter could be slightly 
different in the B° s case and will eventually have to be extracted from experiment. 
It can be extracted easily from the decay channel B® — > DfX. For the decays 
involving a D* we have 



T 



B°DV 



3/- 



(21) 



We shall make the very same assumptions to model Q as we made to model T . 
To do so we have to consider the contribution of b — > ucs to the decay rates of 
the form T(B° -> D S X) which we shall denote by T(B° -> D s X) b ^ uSs . Using the 
same model as the one we used for the wrong charm case in 0, we have 



T(B° S -> D: + X) f 



b—>ucs 
dir ' 



(22) 



where -> D+Xf^'" " 



the current 6 
over a D 



dir — T d ~^. ucs is the direct contribution to this decay over 
ucs. If we take into account the contribution from the decay 
we obtain 



F(B° S -> Dfx) 1 



A "pb — >ucs 
dir 



(23) 



since a _D* + always decays into a We assume that the direct contribution 
T b £ u ~ cs to T(B° S -> D+X) 6 ^ ues is equal to the direct contribution to r(5° -> 
D + X) 6 ^ ucs which is certainly the case in the limit of the heavy quark symmetry. 
We assume that every c quark eventually hadronizes into a D meson. Spin 
counting and isospin symmetry in the decays B — > D*X and B — > Ddi r X through 
the transition b — > ucs allow us to deduce easily that T^ir is equal to 1/8 times the 
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result of the parton calculation. We then obtain the following non-perturbative 
form factors 



QBP.D+ = QBP.D7 = l j 2 (24) 

and 

gB°D:+ =g B°Dr = 3/8 (25) 

In the B° s system, we can use the following approximations 

a = and i^L = J^L = 1. (26) 

M 12 Mi 2 v ; 



These are very good approximations, theory predicts a < 10~ 3 , we then neglect 
the weak phase in p^fj- We see that the only weak phases appearing are e* 7 in 

-,B° B? 



equation (|T^) and e 11 in equation (pTTj) . We can therefore rewrite T Ds s s and 



*° as 



and 



m = n x + e-^n 2 (27) 



B ° = ni + e <7 n 2 , (28) 



where n\ corresponds to the contribution from b —>■ ccs interfering with it-self 
and e* 7 n 2 corresponds to the contribution from b — > cms interfering with b — *> ucs. 
Inserting the formulas for the time dependent decay rates in the definition for 
the CP asymmetry, we obtain 

Acp(t) = -2n 2 sin(AMt)sin( 7 ) 

cn J M 1 (t)2r(5 s ° -> D+X) + M 2 (t)2r(5 s ° -> DjX) + M 3 (t) ' 1 ; 

where Mi(t),M 2 (t) and Ms(t) are given by 

o /AMt\ , o /Am , . 

Ml (t) = cos 2 J + sinh 2 J , (30) 



M 2 (t) = sin 2 (— ) + sinh 2 (— j (31) 



and 



M 3 (t) = -(2nj + 2n 2 ) sinh . (32) 
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The time integrated CP asymmetry is the given by 

— 2n 2 j^2 sin(7) 

AcP = N^yWBQ -+ D+X)+N 2 (x,y)T(B° D~X) + N 3 (x, y) ' (33) 
where Ni(x,y), N 2 (x,y) and N 3 (x,y) are given by 

2 + x 2 v 2 . , 

tfi z,y =T— - + 7^, 34 
1 + x z 4 — y A 

x 2 ?/ 2 

^2 (x, y) = — — j- + - (35) 

and 

N 3 (x, y) = -— y — (2m + 2n 2 cos( 7 )) , (36) 

4 - r 

where we have introduced the parameters x = AM /V and y = AT/T which 
can be measured. At the present time, there is only a lower bound for the 
parameter x in the B° s system namely \x\ > 14 ||. For numerical calculations, 
we shall use \x\ = 20 0. It is also not clear what the actual value of y is, 
for numerical calculations, we shall use the value computed in [||, \y\ = 0.054. 
We set B s = B = 1, m b = 4.8 GeV, m B o = 5.3693 GeV , m c = 1.4 GeV, 
m£) s = 1.9685 GeV, /bo = 210 MeV and Ci = 1. We obtain numerically 

nx = -4.931 %, (37) 
n 2 = -9.526 10" 2 %, 

where rii and n 2 are normalized to the decay width of the B° s meson, using 
r B o = 1.54 10~ 12 s. The decay width needed were computed in ||, we had 

T(B° S -> D+X) = 3.3% (38) 
T(B° S -> £TX) = 64.9%. 

The time integrated CP asymmetry is then given by 

A C p = 1.4 10- 4 sin( 7 ), (39) 

neglecting the term proportional to cos(7) in the denominator. 

3.2 CP asymmetry r(£j(£) vs. r(£J(*) -> D+X) 

In this case the CP asymmetry reads 



^ l J r(fio (t ) _> d 7 x) + r(s°(t) - Dfx) ■ 



(40) 
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B° B° B° B° 

We need to parameterize T t s and T + a . They are given by 

T*l S ° = (T*l B °y = n 1 + e^n 2 (41) 
T Df = { T Dt ) =ni + e^'n 2 . 

The form factors appearing in the transition matrix element were given previously. 
The time dependent CP asymmetry then reads 

A m = -2n 2 sin (AMt) sin( 7 ) 

CP[ 1 M x {t)2Y{B Q s -> D S ~X) + M 2 (t)2T(B° s -> D+Jf) + M 3 (t) ' 1 J 

where Mi(t), M 2 (t) and M 3 (t) were defined in equations ©, © and (||). The 
time integrated CP asymmetry is then given by 

-2n 2T ^sin( 7 ) 



(43) 



VP 



where Nx(x,y), N 2 (x,y) and N s (x,y) are given in equations fl34|), fl35"D and 
The time integrated CP asymmetry is numerically given by 

A C p = 1.4 10- 4 sin( 7 ), (44) 

neglecting the term proportional to cos( 7 ) in the denominator. 

3.3 CP asymmetry T(B° s {t) -> D*JT) vs. r(5°(*) -> L>*X) 

We can without any difficulty deduce from the preceding sections the formulas for 
the CP asymmetries in the one-particle inclusive decays of the form r(B®(t) — > 
D*X) versus T(B°(t) -> L>JX). We have 

m = -3.698 % (45) 
n 2 = -7.145 1(T 2 %, 

which are normalized to the decay width of the B® meson. The decay widths 
needed were computed in 0], we had 

T(B° S -> D* s + X) = 2.5 % (46) 
r(5 s ° -> £>*-X) = 49.6 %. 

The next CP asymmetry we shall consider is defined by 

, m = rcgpg) -» d;-x) - r(gg(t) -» 
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The time integrated CP asymmetry is then given by 



A C p = 1.4 10- 4 sin( 7 ), (48) 

neglecting the term proportional to cos (7) in the denominator. For the asymme- 
try defined by 

, m = rcgpg) -» grx) - r(jg(g -» d;+x) 

^ l J r(s°(t) -> + r(s°(t) -> d^+x) ' 1 J 

we obtain 

^cp = 1-4 10- 4 sin( 7 ), (50) 
neglecting the term proportional to cos(7) in the denominator. 



4 Discussion of the results 

We see that in principle we can extract information on sin (7) from the asym- 
metries calculated in the previous section. But they are small, nevertheless the 
decay widths involved are sizable. It has been proposed to extract information on 
sin(7) from exclusive decays (see e.g. ||) but the decay widths involved are very 
small, typically of the order 10 -4 and one would have to deal with strong phases 
which would make the extraction of sin(7) even more difficult. It could therefore 
be worth to try to extract sin(7) from one-particle inclusive decays in the B® 
system. If the present method is chosen to extract sin (7), it would be interesting 
to test its precision. This could be done by comparing the results obtained for 
sin(2/3) in one-particle inclusive CP asymmetries in the Bd system H with some 
more conventional extraction technique like the "gold-plated" B — > J/^Ks, al- 
though one-particle inclusive CP asymmetries in the B^ system are theoretically 
not as clean as the ones in the B° s system due to the presence of strong phases 

I- 

We still have some large uncertainties, some of them due to the method. The 
corrections to the decay widths could be fairly large, in the worth case of the order 
of 30%. But remember that the decay widths calculated in || are compatible with 
current experimental knowledge. On the other hand, we have large experimental 
uncertainties in the values of x, y and /#o. 

Time dependent CP asymmetries could also allow to extract sin (7), but it is 
not yet clear if the oscillations can be resolved. Predictions depend on the decay 
width difference in the B° s system and it has not yet been possible to measure 
this quantity. 

A way to improve the magnitude of the CP asymmetries would be to do an 
anti-lepton tagging. The decay widths in the denominator are partially respon- 
sible for the low magnitude of the asymmetries, this is particularly true for the 
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CP asymmetries in the B^ system considered in ||. Thus, doing an anti-lepton 
tagging, would slightly improve their magnitude. We would then consider asym- 
metries of the type 

T(B° s (t) -> D S X) NL - T(B° s (t) -> D S X) NL 
^ U T(B»(t)^D s X) NL + T(B*(t)^D s X) NL > { ] 

where NL stands for non-leptonic. In the B® system, we would have time inte- 
grated CP asymmetries of the order 2 10~ 4 sin(7). The factor x = 20 is responsi- 
ble for the low magnitude of the asymmetries. In the system, this effect would 
obviously be larger. 



5 Conclusion 

We have discussed CP asymmetries in one-particle inclusive B° s — > D S X decays. 
The asymmetries are small but would allow to extract sin(7) which is known 
to be difficult. So any new method is probably welcome. It has the advantage, 
in comparison to CP asymmetries of exclusive decays, to have some large de- 
cay widths and in comparison to CP asymmetries of inclusive decays, of being 
experimentally clean. 
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